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ABSTRACT 

We derive constraints that must be satisfied by the sources of ~ 10'^ to ^ 10'** eV cosmic rays, under the 
assumption that the sources are Galactic. We show that while these constraints are not satisfied by ordinary 
supernovae (SNe), which are believed to be the sources of < 10'^ eV cosmic rays, they may be satisfied by the 
recently discovered class of trans-relativistic supernovae (TRSNe), which were observed in association with 
gamma-ray bursts. We define TRSNe as SNe that deposit a large fraction, /« > 10"^, of their kinetic energy 
in mildly relativistic, 7/3 > 1, ejecta. The high velocity ejecta enable particle acceleration to ^ lO'^ eV, and 
the large value of /« (compared to /« ~ 10"^ for ordinary SNe) ensures that if TRSNe produce the observed 
~ 10"* eV cosmic ray flux, they do not overproduce the flux at lower energies. This, combined with the 
estimated rate and energy production of TRSNe, imply that Galactic TRSNe may be the sources of cosmic rays 
with energies up to 10'^ eV. 

Subject headings: acceleration of particles — cosmic rays — gamma-rays: bursts — supernovae: general — 
supernova remnants 
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1. INTRODUCTION 

The spectrum of cosmic rays (CRs) extends from ^ 10^ eV 
to ~ 10^" eV (see lBlandford& Eichleii fT987llAxfordllT994t 
iNagano & Watsonll2000l for reviews). While there is a wide 
consensus that the sources of CRs of energy < 10'^ eV are 
Galactic supernovae (SNe), and that the sources of CRs above 
lO'^ eV are ex tra-Galactic (possibly GRBs or AGN, e.g. 
I WaxmanI l2004d for a recent review), there is no consensus 
regarding the sources of cosmic rays in the energy range of 
-10'^ eVto-lO'^eV. 

The CR spectrum steepens at ~ 10''' eV, the "knee," and 
then extends smoothly to ~ 10"* eV. These characteristics dis- 
favor a transition from Galactic to extra-Galactic sources in 
this energy range, since for such a transition one would nat- 
urally expect either a flattening of the spectrum or a marked 
drop of the flux at the transition energy. In fact, these charac- 
teristics suggest that the same type of sources, or the same 
class of sources, produces the CRs over the entire energy 
range of ~ 10'^ eV to ^ 10'^ eV. Indeed, it has been suggested 
that Galactic SNe are the sources of all CRs in this energy 
range (e.g. Bell & Lucek 2001). 

Over the past several years a new type of supernova explo- 
sions has been discovered, to which we will refer as "trans- 
relativistic supernovae" (TRSNe). We define TRSNe as su- 
pernovae which, unlike ordinary SNe, deposit a significant 
fraction, > 10"^, of their kinetic energy in mildly rela- 
tivistic, 7/3 > 1, ejecta (the corresponding fraction for ordi- 
nary SNe is typically < 10"^; see § l3.2l for a detailed discus- 
sion). Such TRSNe have been discovered recently in associ- 
ation with gamma-ray bursts (GRBs). Hereafter we will use 
the term SNe for ordinary supernovae, which deposit only a 
small fraction of their kinetic energy in relativistic ejecta, and 
the term TRSNe for supernova explosions with /r > 10"^. 

Various terms have been used in the literature to describe 
supernovae associated with GRBs: GRB-SNe, "hypernovae", 
broad-line type Ic SNe, and others. We have chosen here to 
use the term "trans-relativistic supernovae" for several rea- 
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sons. First, it captures the main physical property that charac- 
terizes supernovae that may produce the observed ^ lO'*" eV 
to ^ 10'** eV cosmic-ray flux, namely /r > 10"^. Second, 
supernovae associated with GRBs may not all be similar to 
each other. SN2006aj, for example, is qualitatively differ- 
ent from the other GRB associated supernovae. Although 
there is evidence for this explosion to be characterized by 
/r > 10"^, its total kinetic energy is similar to that of ordi- 
nary SNe, ~ 10''' erg, while the other GRB associated su- 
pernovae appear to be ten times more energetic. Thus, while 
SN2006aj does not fit into the "hypernovae" (i.e. highly ener- 
getic) category, it does fit into the TRSNe category and should 
be considered as a potential source of ^ lO'*" eV to ~ 10'** eV 
cosmic -rays. Finally, it should be kept in mind that not all 
TRSNe are necessarily associated with GRBs. 

We show in this letter that Galactic TRSNe may be the 
sources of - 10'^ eV to - 10'** eV CRs. We first derive in §|2] 
constraints, which apply to both SNe and TRSNe, that must 
be satisfied by Galactic sources of ~ 10'^ eV to ^ 10'** eV. 
We then show in §|3]that while these constraints are not satis- 
fied by SNe, they may be satisfied by TRSNe. Our results are 
summarized and their implications are discussed in § |4] 

2. CONSTRAINTS ON GALACTIC SOURCES OF CRS 

We consider three constraints which must be satisfied by the 
candidate sources of CRs at high energy, e ^ lO'^ eV. First we 
consider in § 12. II constraints on the rate of occurrence of the 
sources in the Galaxy and on their energetics. The constraints 
derived in § 12.11 are applicable to any (Galactic) candidate 
source. Next, we consider in § 12.21 and § 12. 3 1 constraints that 
are applicable to SNe and TRSNe. These sources are expected 
to accelerate particles to high energies through the collision- 
less shocks which they drive into the plasma surrounding the 
exploding star In § 12. 21 we derive constraints on the velocity 
of the shell ejected by the explosion, and in § I2.3l we derive 
constraints on the structure of the ejecta. 

2.1. Rates and energetics 

In order for a certain type of sources to be the main contrib- 
utor to the observed CR flux at some energy eo, the sources 
must satisfy the following conditions: 
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1 . The energy production rate at eq should be sufficient to 
explain the observed flux, and 

2. The Galactic event rate, that is the rate of occurrence 
of the sources in the Galaxy, should not be much lower 
than one event per Galactic confinement time of CRs 
with energy Eq. 

Consider a source with a Galactic occurrence rate of and 
an energy per event of E^. denotes the energy available for 
the acceleration of 10"^ eV CRs, rather than the total energy 
released. In SN explosions, for example, the energy available 
for CR acceleration is the energy deposited in the shock wave 
driven by the ejecta. Since the maximum energy to which par- 
ticles may be accelerated increases with shock velocity (see 
§ 12.2b . the energy available for acceleration of particles to 
10'^ eV is the energy carried by the fastest (outermost) part 
of the ejecta, which drives, at early times, a sufficiently fast 
shock. This energy may be significantly smaller than the to- 
tal kinetic energy of the (slower) ejecta. We assume that the 
source produces an energy spectrum of accelerated particles 
given by 

e^^=ae)E,. (1) 
de 

Here Qi^) is the fraction of Es that is deposited in CRs within 
a logarithmic energy interval around e. The resulting energy 
density of CRs per unit logarithmic particle energy is 



2 ^ , Cs(g)Mre(£)£s 
£ «(£)= 77 : 



(2) 



where n(e) is the number density per unit CR energy, is 
the confinement time of the particles in the Galaxy, and Vc 
is the effective volume of the Galaxy in which the particles 
are confined. We assume that this effective volume is energy 
independent. 

To reduce uncertainties we compare the resulting densities 
with those of low energy, lO' eV, CRs, assuming that in this 
energy range the CRs are accelerated predominantly by SNe. 
Using Eq. (|2|i, the ratio of CR number densities (per unit CR 
energy) at low and high energy, «/,/, = n{e = eij,), is 



A^s Tc(eh) E, Cs(£/<) 
NsN Tc(si) EsN Csn(£/) ' 



(3) 



where issN is the characteristic kinetic energy of a SN explo- 
sion. For £;, = 10"^ eV and £; = 10^ eV the observed ratio is 
(e.g. iNagano & Watsorill2000l and references within) 

^ « (10'Vl0V'(10'Vl0'Y' ° - 10-^ (4) 

efni 

Using Eqs. ^ and (|4]l, we find that the sources of CRs at 
£ = £a = 10''^ eV must satisfy 

-1 
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where A^s =-^s7'c(£a) and A^sn =-^SN'''e(£/) are the average num- 
bers of events that contribute to the observed flux at any given 
time at energies £/, = 10'^ eV and £/ = 10^ eV respectively, 
and C = CsClO'** eV)/CsN(10'' eV). Using the confinement 



time of 1 GeV protons in the Galaxy, tc(W e V) « 10'^ yr 
(IWebber & SoutoulllToOSl: lYanasak et al.ll200l . Eq. (01 may 
be written as 



.0.03C-'(|^ 



N, 



SN.-2, 



(6) 



where Nsn = O.QINsn-2 yr '- The requirement that the event 
rate should exceed one event per Galactic confinement time, 
Na^l, implies 

^ A < 0mNsN,-2- (7) 

The confinement time of CRs of energy ~ lO'^ eV is likely 
to be significantly larger than the light crossing time of the 
Galaxy, ^ 10"*yr, due to the following argument. The ob- 
served energy evolution of the depth of maximum, X^ax, 
of extreme air showers suggest that composition of CRs at 
~ lO^ ^eV is dominated by nuclei with moderate at omic num- 
ber ZljOaiiiEilillillllNaiinE&Mio^^ Consider 
therefore the Larmor radius of a CR particle with atomic num- 
ber Z = lOZi and an energy e = 10'^£i8eV, 

Rl ~ 40B-_lsSi&Zi' pc, 

where B = 3B-5.5/iG is the Galactic magnetic field. Since 
the Larmor radius of such particles is much smaller than 
the thickness of the CR disk, ^ 1 kpc, the propagation of 
£ = lO'^eV CRs is expected to be strongly affected by the 
Galactic magnetic field. The transition to non confined par- 
ticles, for which the confinement time is comparable to the 
light crossing time, should then occur at > 10 Zi eV. Using 
therefore rc(10'^ eV) = r4,5 lO'^-^yr, Eqs. Q and @ give 



A^s 



IO-V7' 



— ) r'. 

Esn' ^ 



(8) 



The constraint that the Galactic event rate be larger than one 
event per confinement time can be written as: 



A^s 



> 3 X lO-'u',Nsk-2- 



(9) 



Eqs. ([7]), ([8]l and (|9]l (any one of these can be derived from the 
other two) describe the constraints on the occurrence rate and 
energetics of the sources of ^ lO'^ eV CRs. 

A note is in order regarding the value of the parameter 
C = C(10'^ eV)/CsN(10'') eV. In the following, we consider 
sources that accelerate particles through collisionless shocks 
in the surrounding medium driven by fast ejecta. Note that al- 
though a complete theory, based on first principles, of particle 
acceleration in collisionless shocks is not available, there are 
strong theoretical arguments and observational evidence indi- 
cating that for both non-relativistic and relativistic shocks a 
considerable fraction of the post shock energy can be con- 
verted to relativistic par ticles with a spectrum that follows 
e^dN/de w £" (see, e.g. iBlandford & Eichleiill987l: lAxfordl 
iT994; Waxman 2006| for reviews). If both SNe and the 
sources of 10'^ eV CRs deposit a significant fraction of their 
energy in accelerated particles with a spectrum dN /de oc £~^, 
the value of Q would be of order unity. 

2.2. Maximum energy 

We assume that the particles are accelerated in the vicinity 
of a collisionless shock wave with (a time dependent) veloc- 
ity Vs = j3sc that is driven into the surrounding medium by an 
ejecta of mass Mgj, initial velocity Vgj = /3ejC and kinetic en- 
ergy ^ej = (7ej - l)MejC^. Here = (1 - /Jgj)'/^ is the ejecta's 
Lorenz factor We consider two options for the density dis- 
tribution of the surrounding medium: A homogeneous inter- 
stellar medium (ISM) or a progenitor wind with a density pro- 
file p oc r~^. 



3 



At any given time, the maximum energy to which particles 
can be accelerated is given by 



' ZePsBdrs 



(10) 



(e.g. IWaxmanI I2OO4H for a recent review). Here, rs is the 
radius of the shock and is the minimum of the mag- 
netic field in the pre-shock (upstream) and post-shock (down- 
stream) plasma, measured in downstream frame. We should 
stress that the energy given by Eq. (fTol l is an upper limit, and 
that the maximum energies attained are likely to be lower. 
For example, if particles are accelerated via Fermi's diffusive 
shock acceleration, the maximal energy would be a few times 
lower for Bohm limit diffusion, and much lower for slower 
diffusion. 

There is growing evidence that the magnetic field in the 
vicinity of non-relativistic collisionless shocks in SNRs and 
gamma-ray bursts (GRBs) afterglows is amplified to values 
greatly exceeding the ambient magnetic field both in the up- 
stream and downstream. In the past few years, high resolu- 
tion X-ray observations have provided indirect evidence for 
magnetic fields of values as high as 100 /iG in the (vs ^_ 
few X 1000 km s~') shocks of young SNRs (s ee lBamba et al.l 
I2003t IVink & La ming 2003t IVolk et ani2005 >). For relativis- 
tic shocks in the afterglow phase of GRBs, a very large mag- 
netic fi^ldJnthe_downstream is inferred from measurements 
(see IWaxmanll200 6*. for review), that can not be explained as 
a mere compression of the ambient field and requires a fair 
fraction of the kinetic energy of the shock to be transformed 
into magnetic energy. Li & Waxman (2006) have furthermore 
shown that in GRB afterglow shocks strong amplifications of 
the magnetic field in the upstream can be also inferred. We 
therefore assume that the magnetic field is amplified by the 
shock both in the upstream and in the downstream to values 
close to equipartition 



eB87r/9Vj7| 



(11) 



Initially, the shock and ejecta velocities (which are similar) 
are constant with time while the radius increases linearly. For 
the density distributions considered here the maximal accel- 
eration energy is attained when the ejecta begins to deceler- 
ate. This happens roughly when the energy in the swept up 
material is comparable to the initial ejecta energy, i.e. when 



swept 



' Mej/7ej. In case the shock is propagating into a ho- 



mogeneous ISM, the deceleration radius is given by 



1/3 



rho 



(12) 



where no is the ISM particle density per cm^. Using Eqs. ( fTOl i. 
(fTTI ) and ( fT2] l. the maximum energy of an accelerated particle 
is 



£<3xlO'^Zie^^i 



V lOMp 



1/3 



.7feV. 



(13) 



Here = O.leg-i and /3ej = 10 ^/3ej -2- 

For a shock propagating into a progenitor wind with a den- 
sity profile p = M/47rr^v,v, where M is the mass loss rate and 
Vh, is the wind velocity, the deceleration radius is 



^ wind ' 



7ejM ' 



(14) 



Using Eqs. ( fTOl i. (fTTT l and ( fT4] l. the maximum energy of an 
accelerated particle is 



e< lO'^Zie^/?, 



M_5 



1/2 



/3ej.-27ejeV. 



(15) 



Here M = IO'^M-sMq/yt and v^^. = 10'^v„,.8cm/s. The maxi- 
mum energy is independent in this case of the ejecta mass, i.e. 
of the deceleration radius, since for a 1 /r^ density profile the 
maximum energy (for a fixed shock velocity) is independent 
of the shock radius. 

Note, that this estimate is valid provided shock deceleration 
occurs while it is propagating within the p oc wind profile. 
The p oc dependence of wind density extends up to a radius 
where p,^ vf^, ~ Pism, ™d the wind mass contained within this 
radius is 



(4^v„P„m)'/2 



; 0.3M. 



l/2pl/2 ■ 



(16) 



where the ISM pressure is Pjsm = 0.3eV/cm^f/sM.-o.5- The 
maximum energy is given by eq. ( fTSI ) for Mej/jej < M^,, and 
it is lower otherwise. 

2.3. Energy distribution within the ejecta 

SN explosions produce non-uniform ejecta, with velocity 
rising towards the front edge of the ejecta. The shock driven 
by the fastest, outermost part of the ejecta is capable of ac- 
celerating particles to the highest energy. Here we consider 
the case where a small and fast part of the ejecta is respon- 
sible for accelerating cosmic rays to energies £/, ^ 10'^ eV, 
and constrain the energy distribution of the slower parts of the 
ejecta by requiring that they do not produce a flux of lower 
energy CRs that exceeds the observed one. 

We denote the amount of energy in ejecta moving faster 
than (3c by Ek{> (3). Since we consider in this section a lim- 
ited range of particle energies, we assume that each part of 
the ejecta contributes a constant fraction ^ej to each logarith- 
mic bin of particle energies at e < emax(/3) (ignoring the pos- 
sible energy dependence of ^ej)- Under this assumption, the 
resulting energy flux per logarithmic particle energy interval 
is eMe) oc Tc{e)Ek{> /3(e)). Here (3(e) (oc e'/^ foj. < i) 
is the shock velocity required for accelerating up to e, i.e. 
emax(/3) = e. 

Assuming that the part of the ejecta with velocities larger 
than some (3h is responsible for the CR spectrum at some e/,, 
requiring the slower parts of the ejecta do not produce a flux 
of cosmic rays at lower energy that exceeds the observed flux 
implies 



Ek{>(3{e)) ^ e2«obs(£) 



Ek{> (3{eh)) elnobsish) 



TciSh) 



(17) 



where «obs stands for the observed CR number density (per 
unit CR energy). Since the measured energy flux per log- 
arithmic particle energy interval in the range lO'^eV < e < 
10'^ eV is roughly proportional to E~\ and the confinement 
time is expected to decrease with energy, the energy distribu- 
tion must satisfy 



EkO (3(e)) 
EkO /?(£/,)) 



< 



(18) 
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3. APPLICATION TO SNE AND TRSNE 
3.1. SNe 

Eqs. (fTTt and ( fTSl ) imply that typical SNe driven shocks, 
characterized by (3ej ~ 10"^, are probably too slow to acceler- 
ate particles to 10'** eV, especially considering that these are 
optimistic upper limits. Another difficulty in attributing the 
production of 10'** eV CRs to ordinary SNe is related to the 
constraint of Eq. (O. According to this constraint, if SNe 
were capable of accelerating particles to energies as high as 
10'^ eV, the flux produced at these energies would exceed the 
observed flux by a large factor unless C ^ lO""*, i.e. unless the 
energy carried by 10"* eV CRs accelerated by the SN shock is 
smaller by a factor of 10'* than that carried by 10^ eV CRs ac- 
celerated by the same shock. Such a strong suppression is not 
expected for the case when the maximum acceleration energy 
exceeds lO'** eV, since in this case we expect equal energy to 
be deposited per logarithmic CR energy interval (see discus- 
sion at the end of section 32. 11 1. For example, attributing this 
factor to an acceleration spectrum dn /c/e cx £~p with p differ- 
ent from 2, would require /:> > 2.4 which seems inconsistent 
with theory and SNR radio observations. 

One may argue that the above problem may be circum- 
vented by assuming that lO'^ eV CRs are produced only by 
the fast, [3 3> 10"^, component of the SN ejecta, which car- 
ries only a small fraction of the total kinetic energy. This 
would lead, however, to violation of the constraint of Eq. (fTsT l. 
The typical energy distribution in SN ejecta is Ek{> f3) oc /3"^ 
for BSG progenitors and Ek( > (3) oc /3"^ for RSG progeni- 
tors (iMatzner & McKeelll999h . At non relativistic velocities 
Emax OC 0^ , implying /3(e)) oc e~^l^ or Ek{> /3(e)) oc 
for BSGs and RGBs respectively. This is in clear contradic- 
tion with Eq. (fTsT l. That is, if one assumes that the fast part 
of the ejecta produces the observed 10''^ eV CR flux, then the 
predicted flux at lower energies would far exceed the observed 
flujL 

Bel l & Lucekl (1200 1) have argued that SN driven shocks are 
candidate sources of CRs above lO'** eV, based on an esti- 
mate of the maximum acceleration energy. They have con- 
sidered a shock of velocity of 40,000km/s driven into a mas- 
sive wind, characterized by Ms ^ 1 and v,v 8 ^ 10"^, as in- 
ferred from radio observations of SN1993J. Indeed, for such 
parameters the SN driven shock may allow acceleration to 
> 10'^ eV, as indicated by Eq. ( fTsT l. However, it should be 
realized that SN1993J is still young, and that the measured 
velocity therefore represents only the fast front edge of the 
ejecta, which contains only a small fraction of the ejecta mass 
an d energy. This i s well illustrated in the modeling of this SN 
by IWoosley et al.l (1 19941) . where only a small fraction of the 
total mass and energy lies in such a high velocity component 
of the ejecta (see e.g. their fig. 6). As explained above, if 
the tail of the ejecta velocity distribution can account for the 
cosmic ray flux at some energy e;, ^ 10"* eV, then the steep 
profile of energy as a function of velocity would imply a CR 
flux at lower energies that far exceeds the observed flux. 

3.2. TRSNe 

During the past decade, four GRBs were observed 
to be associated with core co llapse SNe o f typ e 
Ic: GRB980425/SN199 8bw (iGalama et all 119981) 
GRB030329/SN2003dh (Hio rth et al.1 l2003t IStanek et al.l 
120031) . GRB031203/SN20031W (Tagliafem et al.' "200?), 
and XRF060218/SN2006aj (Campana et al. 2006) . All but 
GRB030329 were very faint, 4-5 orders of magnitude fainter 



than an average cosmological GRB. In all 4 cases there is 
evidence from radio observations for deposition of a signifi- 
cant part of the kinetic energy in a mildly relativistic part of 
the ejecta (ISoderberg et al.ll2006l and references therein). Al- 
though radio observations alone do not allow one to uniquely 
determine the energy carried by the relativistic ejecta 
dWaxman & Loeblll999t rWaxmanll2004ah . additional obser- 
vational information enabled such a unique determination 
in two cases. For GRB980425/SN1998bw there is a robust 
estimate of the energy deposited in the fast, /3 > 0.85, part of 
the ejecta, based on long term radio and X-ray observations, 
which yield £(/3 > 0.85) = lO'*''-^ erg (Waxman 2004b). 
Similar values are infeiTed for XRF060218/ SN2006ai from 
prompt and afterglow X-ray observations (IWaxman et al.l 
120071) . 

The SN fight curves of SN1998bw, SN2003dh, and 
SN20031w appear to be qualitatively similar, implying SN 
ejecta mass and energy of Egj ^ 5 x 10^^ erg, Mej ^ IOM0 
(e.g. table 6 of iMazzaU et al.ll2006ah . Note, that the energy 
may be lower, ^ 1 x 10^^ erg, if the explosion is not spheri- 
cal, which is likely (iMaeda et a l. 2006). The SN light curve 
of SN2006aj is qualitatively different than the others, imply- 
ing Eej ~ 2 X 10^' erg and Mej — 2Mq (in thi s case, the ex- 
plosion is probably not highly non spherical; iMazzali et al.l 
I2006bl) . The infeiTed E,j and E(f3 > 0.85) imply that a sig- 
nificant fraction of the kinetic energy is deposited in a fast, 
mildly relativistic part of the ejecta, £'(/3 > 0.85) /E^j > 10"^. 

Two points should be emphasized here. First, GRBs associ- 
ated with TRSNe are different than cosmological GRBs: The 
amount of energy carried by relativistic ejecta in TRSN explo- 
sions, ~ lO'*'''^ erg, is much lower than that of cosmological 
GRBs, which is 10^' erg, and the fastest part of the ejecta 
in a TRSN is only mildly relativistic, unlike the highly rela- 
tivistic ejecta inferred from observations of GRBs. Second, 
the mechanism responsible for depositing ^ lO'*'''^ erg in the 
mildly relativistic part of the ejecta is not understood. For the 
inferred values of Eej/Mejc^ ^ 10"^, the acceleration of the 
SN shock near the edge of the star is expected to deposit only 
^ lQ~^Eej in th e part of the ejecta expanding with 7/? > 1 
("Tan e t all 1200 Ih . This suggests that the mildly relativistic 
component is driven not (only) by the spherical SN shock 
propagating through the envelope, but possibly by a more rel- 
ativistic component of the the explosion, e.g. a relativistic 
Jet pro pagating through the star (lAloy et al.l2000l : IZhang et all 
l2003h . For the present analysis, we use therefore only the ob- 
servational constraints on the energy distribution within the 
ejecta. 

With only 4 detected events, the TRSNe rate is rather uncer- 
tain. Estimates of the local (z = 0) rate range from ~ 10^ to ~ 
lO^Gp c"^ yr"' ([G uetta & Delia Valle' "2007; Soderberg^taD 
2006; Eiang et a l. 2006; Pian et al. 2006). Comparing this 
to flie local SN rate, w lO^Gpc"^ yr"' (e.g. ICappellaro et alJ 
119991) . we find A'trsn/A'sn » 10"^-^^° ^ This estimate is con- 
sistent with the estimate of iPodsiadlowski eT al. ( 2004) of the 
rate of "hypernovae" in galaxies similar to the Milky Way. As 
mentioned in the introduction, the TRSNe rate may be higher 
if some TRSNe are not associated with GRBs (and hence were 
not identified following a GRB trigger). With these estimates 
at hand, we are now ready to determine whether or not TRSNe 
meet the constraints derived in § |2] 

First, for /3 w 0.8 Eqs. ( fT3] l and dTsl l imply that the mildly 
relativistic part of TRSNe ejecta may indeed accelerate CRs 
up to energy exceeding 10'^ eV. Second, for an energy 
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Er ^ 10^^^ erg in the relativistic part of the ejecta we have 
iis/^SN ^ 10"' satisfying the energetics constraint, Eq. (|7]l, 
and ^trsn/^sn ~ io~^'^*°'^ satisfying the rate constraint, 
Eq. ([8j. Note, that for TRSNe we expect ( ^ I, based on their 
radio and X-ray observations. These observations imply that 
the coUisionless shocks driven by the mildly relativistic ejecta 
into the wind surrounding the progenitor transfer a significant 
fraction of the energy to a population of shock accelerated 
electrons, and that the accelerated electron spectrum follows 
e^dN/de w e° over a wide range of energies, from electrons 
emitting synchrotron radiation in the radio band over at least 
~ 5 deca des of energy to t hose emitting radiation in the X- 
ray band (' Waxmanl l2004b !) . Thus, if nuclei are accelerated 
with similar efficiency to a similar energy power spectrum, the 
fraction of energy deposited in ~ 10''^ eV CRs would be sim- 
ilar to that deposited by SNe in CRs accelerated to 10^ eV. 

Let us next consider the constraint of Eq. ST% on the energy 
distribution within the ejecta. Consider a TRSN with a total 
kinetic energy of Efj = lO^^Egj 52 erg, an ejecta mass Mej = 
lOMej.iM0 and an energy Er in a mildly relativistic, f3 > 0.8, 
part of the ejecta. Using Eq. (fTsT l. the ratio of the maximum 
energies of particles accelerated by the relativistic part of the 
ejecta, £mnx,r, and by the bulk of the ejecta, emax.c/, is 

emax.,-/emax,<.i « 10^''£',]'52Mej,i , (19) 

and the constraint of Eq. ([TtT i is satisfied in the energy range 
of 10'^ eV-lOi'^eV as long as 



7^ i/-i49^ej.52 Tci^m-dx,!-) 

' M '( 



Assuming that Tc oc e we have 



£^>1049^103^erg. 



(20) 



(21) 



ej.l 



An upper limit to the value of S in the relevant energy range, 
lO'^ eV- 10'^ eV, may be obtained as follows. First note that 
the confinement time of Z = 10 particles at e = lO^ZeV is 
10^-^ yr while and at e = lO'^Z eV it is larger than 10^^ yr, 
implying that the average value of S between these energies 
satisfies 

log(e/,/e/) 



S = 



<0.4. 



log[Tc(eh)/Tc(ei)] 

In fact, the value of 5 in the energy range of 10'^ eV- 10'^ eV 
must be lower than this average. The grammage (column 
density) traversed by CRs of energy < 1 TeV before they es- 
cape o ur Galaxy is Sconf ~ 9 ( £'/10Z GeV)~'^g cm~^ with S ~ 
0.6 (Engelmann et alj [1990t [Stephens & Streitmatteri [19981: 
[Webber et al. 2003), suggesting cx 10^-^(£/Z GeV)"" '^ yr 
for E/Z < 10'2 eV. As the value of djit energies lO^Z eV- 
lO'^Z eV is higher than the average, S < 0.4, the value of S 
must be lower at energies lO'-'Z eV- lO'^Z eV. It is reason- 
able to assume therefore lO^*"^ < 10, which implies that the 
constraint of Eq. (ISTT i may be satisfied for reasonable values 
of the parameters. 

Finally, we note that i Milgrom & Usovl (Il996h have sug- 
gested that cosmological GRBs, which occur in our galaxy 
once every ~ 10^ years, could be the sources of the CRs 
in the energy range of 10'"* eV to 10'^ eV. This sugges- 
tion is qui te distinct from the one presented in this arti- 
cle. While iMilgrom & Usovl (11996!) considered acceleration 
by the highly relativistic jets of cosmological GRBs, where 

10''' erg is carried by a highly relativistic, 7 10^'^ ejecta, 
we consider acceleration by mildly relativistic, 7/3 ~ 1, less 
energetic, ^ 10"*^^ erg, and more frequent TRSN explosions. 



4. DISCUSSION 

We derived constraints that must be satisfied by the sources 
of ~ 10'^ to 10'^ eV cosmic rays, under the assumption that 
the sources are Galactic and that the CRs in this energy range 
are confined by the Galactic magnetic field (as explained in 
§ 12.11 the CRs in this energy range are likely to be confined 
by the Galactic magnetic field since their Larmor radius is 
much smaller than the CR disk thickness). In § |2.1| we have 
shown that the Galactic occurrence rate, .^Vs, of sources pro- 
ducing ^ 10'^ eV CRs and their energy production E^ per 
event must satisfy the constraints given by Eqs. d?) and dHJ 
[Efi is the energy available for the acceleration of 10'^ eV 
CRs, rather than the total event energy, as explained preced- 
ing eq. ([U]. In particular, E^ must be considerably lower than 
the kinetic energy of a typical SNe, E^ < 0.03Nsn-2C~^ Esf^, 

where Nsj^ = 0.01Nsn.-2 yr~' is the Galactic SNe rate. ( in 
these equations is the ratio of the fraction of £5 deposited in 
10'^ eV CRs to the fraction of the kinetic SN energy that is 
deposited in ^ 10'' eV CRs. These constrai nts m ust besatis- 
fied by any candidate Galactic source. In § 12.21 and § I2.3l we 
derived additional constraints that are applicable to SNe and 
TRSNe. These sources are expected to accelerate particles to 
high energies through the coUisionless shocks that th ey dr ive 
into the plasma surrounding the exploding star In § I2.2l we 
derived constraints on the velocity of the shells ejected by the 
explosions [eqs. (fTsT i & (fTsTll, and in § 12.31 we derived con- 
straints on the structure of the ejecta [eqs. ST7\ & (fTSTll. 

In § 13.11 we have shown that ordinary SNe are unlikely to 
be the sources of 10'^ eV CRs, since they are unlikely to be 
able to accelerate particles to such energy and if they could, 
they would produce fluxes far exceeding the observed flux at 
these energies. We have furthermore shown that assuming 
that the flux of 10"^ eV CRs is produced by the fastest part 
of the ejecta of SNe, which carry a small fraction of the ejecta 
energy, would predict a CR flux at lower energy, ~ 10'^ eV, 
which far e xceeds the observed flux. 

In § 13.21 we have shown that Galactic TRSNe may be the 
sources of high energy, ^ 10'^ eV, CRs. The mildly relativis- 
tic shocks driven by such ejecta are likely to be able to accel- 
erate particles to energies exceeding lO'^ eV. The estimated 
rates of TRSNe combined with the typical energy they deposit 
in mildly relativistic ejecta yield a flux of 10''^ eV CRs which 
is comparable to the observed flux (under the assumption that 
the efficiency of the acceleration of particles in these shocks 
is similar to that of SNe and that the accelerated particle spec- 
trum is close to dN/de oc e~^, which, as explained in § l3.2l is 
supported by radio and X-ray observations). The large frac- 
tion, > 10"^, of the kinetic energy deposited by such explo- 
sions in mildly relativistic ejecta ensures that if TRSNe are 
responsible for producing the observed flux at lO'^ eV they 
do not overproduce lower energy CRs. Furthermore, since 
the TRSNe are a subset of the larger SNe population, their 
CR production is expected to smoothly join the lower energy 
CRs produced by SNe. 

This research was partially supported by AEC, Minerva and 
ISF grants (RB, BK, EW). EW and BK thank the Institute for 
Advanced Study for its hospitality during a period where part 
of this research was carried out. AM acknowledges support 
from the Frank and Peggy Taplin membership at the Institute 
for Advanced Study. While this paper was in preparation, we 
became aware of work in preparation by X. Wang, S. Raz- 
zaque, R Meszaros and Z. Dai proposing hypemovae as cos- 
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mological sources of high energy CRs dWang et alj|2007h . In the present work we have focused on Galactic sources. 
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